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We consider the particular class of f{R) gravities minimally coupled with Yang - Mills (YM) field 
in which the Ricci scalar — Ro = constant in all dimensions d > 4. Even in this restricted class 
the spacetime has unlimited scopes determined by an equation of state of the form Peff = <^P- 
Depending on the distance from the origin (or horizon of a black hole) the state function uj (r) takes 
different values. It is observed that — > i (the ultra relativistic case in 4 - dimensions) and a; — > — 1 
(the cosmological constant) are the limiting values of our state function cj (r) in a spacetime centered 
by a black hole. This suggests that having a constant to throughout spacetime around a charged 
black hole in f{R) gravity with constant scalar curvature is a myth. 
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^ • I. INTRODUCTION 

For a number of reasons, ranging from dark energy and accelerated expansion of the universe to astronomical 
. tests, modified version of general relativity gained considerable interest in recent times. f{R) gravity, in particular, 
attracted much attention in this context (see [l| for comprehensive reviews of the subject). The reason for this trend 
may be attributed to the dependence of its Lagrangian on the Ricci scalar alone, so that it can be handled relatively 
simpler in comparison with the higher order curvature invariants. Depending on the structure of the function f{R) the 
nonlinearity creates curvature sources which may be interpreted as 'sources without sources', manifesting themselves 
I ■ in the Einstein equations. Beside these curvature (or geometrical) sources there may be true physical sources that 
^jjl contribute together with the former to determine the total source in the problem. It should be added that owing to 
I— I. highly nonlinear structure of the underlying field equations attaining exact solutions is not an easy task at all. In spite 
of all odds many exact solutions have been obtained from ab initio assumed f{R) functions. To recall an example 
, we refer to the choice f{R) = R^ , {N = an arbitrary number) which attains an electromagnetic - like curvature 
• source, so that N ^ 1 can be interpreted as an 'electric charge without charge' [11 . That is, the resulting geometry 
(y-^ [ becomes equivalent to the Reissner-Nordstrom (RN) geometry in a spherically symmetry metric ansatz of Einstein's 

■ gravity. This particular example reveals that the failure of certain tests related to Solar System / Cosmology in f{R) 
[ gravity is accountable by the curvature sources in the Einstein Hilbert action. Equivalence with f{R) = i?+ (scalar 
. fields) provides another such example beside the electromagnetic one. More recently we obtained a large class of 

■ non-analytical f{R) gravity solutions minimally coupled with Yang- Mills (YM) field [3|. Even more to this the YM 
T-H , field was allowed to be a nonlinear theory in which the power- YM constitutes a particular example in all higher 
^-H ■ dimensions. In particular, in d = 6, f{R) = VR solves the Einstein- Yang-Mills (EYM) system exactly. For d = 4 

^ . our solution for nonabelian gauge reduces to an abelian one which may be considered as an Einstein-Maxwell (EM) 
• 'H ■ solution [1]. 

^\ I Previously f{R) gravity coupled non-minimally with Yang-Mills and Maxwell matter sources have been considered 
0. In this paper we consider a particular class within minimally coupled YM field in /(i?) gravity with the conditions 
that the scalar curvature R = Rq = constant and the trace of the YM energy-momentum tensor is zero. (To see 
other black hole solutions with matter in f{R) gravity we refer to Ref. Q). Contrary to our expectations this 
turns out to be a non-trivial class with far-reaching consequences. Our spacetime is chosen spherically symmetric 
to be in accord with the spherically symmetric Wu-Yang ansatz for the YM field. The field equations admit exact 
solutions in all dimensions d > 4 with the physical parameters; mass (m) of the black hole, YM charge {Q) and the 
scalar curvature {Ro) of the space time. In this picture we note that cosmological constant arises automatically as 
proportional to Rq. From physics stand point, considering equation of state in the form Peff = ^jjp, with effective 
pressure (Peff) and density (p), important results are obtained as follow: For a critical value oi r = we have 
— 1 < w (r) < for r > Tc and < w (r) < ^ for r < Vc- Remarkably this amounts to a sign shift in the effective 
pressure to account for the accelerated expansion of a universe centered by a charged black hole. In general the critical 
distance is thermodynamically unstable so that the universe undergoes the phase of accelerated expansion beyond 
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that particular distance. Absence of the phantom era (i.e. uj < —1) is also manifest. Alternatively, in the limit r — > 
it yields a; — > i (i.e. 4— dimensional ultra relativistic case), while for r — > oo we have uj — > —1, the case of a pure 
cosmological constant. Let us note that the latter case corresponds to vanishing of the YM field. Stated otherwise, in 
the overall space time we do not have a fixed value for w. Depending on the distance from the center (or horizon) of a 
black hole we have a varying state parameter ui (r). The same argument in the Friedmann- Robertson- Walker (FRW) 
version of the theory implies that beyond a critical time t = tc, oj (t) changes its role and a different type of matter 
becomes active. It is known that for the normal and dark matters which provide clustering both the weak energy 
condition (WEC) and the strong energy condition (SEC) must be satisfied. In the case of dark energy on the other 
hand WEC is satisfied while SEC is violated. In Appendix A we analyze the energy conditions thoroughly covering all 
dimensions. Although our metric ansatz is chosen to be spherically symmetric so that the constant scalar curvature 
J?o > 0, in order to prepare ground for the topological black holes we consider the case of i?o < as well. 

Organization of the paper is as follows. In Sec. II we introduce our formalism and present exact solutions. The 
analysis of our solution with thermodynamical functions is considered in Sec. III. We complete the paper with 
Conclusion which appears in Sec. IV. 



II. THE FORMALISM AND SOLUTION FOR R =CONSTANT. 



We choose the action as (Our unit convention is chosen such that c = G = 1 so that k = Stt) 

7(«) 



/A, 



2k 



+ £{F) 



(1) 



in which / (i?) is a real function of Ricci scalar R and L (F) is the nonlinear YM Lagrangian with F = 
. Obviously the particular choice £ {F) = —-^F will reduce to the case of standard YM theory. The 



YM field 2— form components are given by 



2 '"^ 



(2) 



with the internal index (a) running over the degrees of freedom of the nonabelian YM gauge field. Variation of the 
action with respect to the metric gives the EYM field equations as 



in which 



= £ (F) 6; - tr j (F) , 

dC{F) 



(3) 

(4) 



Cf{F) 



dF 



Our notation here is as follows: fa = D/ij = ^ix^^fn = —^9^ {\/—gd'^) /r , R'j^ is the Ricci tensor and 



The trace of the EYM equation (3) yields 



pm 



fRR + {d-i)nfR--f = KT 



in which T = TI^. The SO {d — 1) gauge group YM potentials are given by 

^(a) ^ .^cgj^.^ x'dx^, Q = YM magnetic charge, = J^^h 



d-l 



(5) 
(6) 

(7) 



i=l 



2 < i + 1 < i < d - 1, and 1 < a < (d - 2) (rf - 1) /2, 
Xi = rcos6d-ssm6d-4...sm6i, 0:2 = rsin0(i_3sin^(i_4...sin0i, 
xs = rcos^d_4sin^d_5...sin^i, X4 = rsin^(i_4sin^(i_5...sin^i. 



Xd-2 = r cos 61, 



in which cj^},^^ are the non-zero structure constants of ('^ ^K*^ ^) —parameter Lie group 5 0,11] ■ The metric ansatz 



is sphericahy symmetric which reads 



ds' = -A (r) dt' + -—+ r'dnl_^, 
A{r) 



with the only unknown function A (r) and the sohd angle element 



dni_2 = del + EU sin" del 

i=2j = l 



with 



< ed-2 < 27r, < 6*, < TT, 1 < i < d - 3. 
Variation of the action with respect to A^"' implies the YM equations 

d *F^'''>Lf (F) + ^cI^J(^^^Lf (F) aC") a* F^'^) = 0, 
in which cr is a coupling constant and * means duality. One may show that the YM invariant satisfies 



F 



-tr ( F^°-' 



fa)p(aW^ ^ {d-2){d-3)Q^ 
4r4 



and 



while 



tr 



(d-3)( 



which leads us to the exact form of the energy momentum tensor 



= diag 



(d-3)( 



-Cf,C- 



{d - 3) 



{d~3)Q^' 



Here the trace of Tj^ becomes 



T = T'f = dC- AFCf 



and therefore with Eq. (3) we find 



/ = ^ UbR + (d - 1) DfR ~n{dC- AFCf)] ■ 
To proceed further we set the trace of energy momentum tensor to be zero i.e., 

dC-AFCp = 

which leads to a power Maxwell Lagrangian 

1 d 
C ^ Fi. 
47r 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



Here for our convenience the integration constant is set to be On the other hand, the constant curvature R — R^, 
and the zero trace condition together imply 



/' (i?o) Ro - 2/ (^0) = 0. 



(19) 
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This equation admits 

f{Ro)^Rl, (20) 

where the integration constant is set to be one. One can easily write the Einstein equations as 

g; = ^ f; (21) 

where 

^ " /(i?0)d K ^ > 

(d-2)i?o 

"^^^^ - 2d~' ^^^^ 

and in which is given by (4). The constancy of the Ricci scalar amounts to 

_ r^A" + 2{d-2)rA' + {d-2) [d-i) {A - 1) ^ 

which yields 

^^^^ d{d-\f^ ^^^^ 
where a and m are two integration constants. From the Einstein equations one identifies the constant a as 

d{d-2)R~ V 4 y 

In the next section we investigate physical properties of our solution in all dimensions. 

III. ANALYSIS OF THE SOLUTION 
A. 4— dimensions 

1. Thermodynamics 

In 4— dimensions, we know that the nonabelian 50(3) gauge field coincides with the abelian U{1) Maxwell field Q]. 
Due to its importance we shall study the 4— dimensional case separately and give the results explicitly. First of all, 
in 4— dimensions the metric function becomes 

and the form of action reads as 



S — I d Xy/—g 

in which 



2k 



C{F) 



(28) 



f{R)^R\ (29) 

with R = Rq and 

^iF)--^F- (30) 



By assumption, Rq gets positive / negative values and the resulting spacetime becomes de-Sitter / anti de-Sitter, 

Ro 

4 



type in f{R) = R^ theory respectively, with effective cosmological constant A^fj = Let us add that in order 
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to preserve the sign of the charge term in (27) we must abide by the choice Rq > 0. However, simultaneous hmits 
— and Rq — >■ 0, so that ^ = Ao =constant, leads also to an acceptable solution within / (R) gravity 0]- It is 
not difficult to see here that m is the ADM mass of the resulting black hole. Viability of the pure f(R) = R^ model 
which has recently been considered critically 10] is known to avoid the Dolgov-Kawasaki instability Further, in 
the late time behaviour of the expanding universe (i.e. for r — >■ oo) it asymptotes to the de Sitter / anti de Sitter 
form. With reference to [13] we admit that sourceless /(i?) = R^ model doesn't possess a good record as far as the 
Solar System tests are concerned. Herein we have sources and wish to address the universe at large. Now, we follow 
[l^ to give the form of the entropy akin to the possible black hole solution. From the area relation the entropy of the 
modified gravity with constant curvature is given by 



which upon insertion from (19) becomes 



S = ^f'iRo) (31) 



^=^/(i?o) = 27ri?or^ (32) 



where indicates the event horizon. The Hawking temperature and heat capacity are given respectively by 



and 



r -T 4^-^o^^ {^14 ^Rprl + 2Q^) 

^"dTn (Rlrl + ARorl-GQ^) ' ^^^^ 

Here we note that for the case of zero YM charge {Q — 0) one finds 

^'3 = ^-9I^^'*'^(4TlvD ^''^ 

which clearly shows from (34) that for i?o > 0, the YM source brings in the possibility of having a phase change. This 
is depicted in Fig. 1. In Fig.s lA and IC we plot the horizon radius versus mass for Q — Q and Q — 1. Similarly in 
Fig.s IB and ID we plot the heat capacity C for Q = and Cq for Q = 1 to see the drastic difference. It is observed 
that for (5 = (Fig. IB) the heat capacity is regular whereas for Q = 1 (Fig. ID), Cq is a discontinuous function 
signalling a phase change. 

2. Energy Conditions 

From the energy conditions (see Appendix, A2) the density and principal pressures are given as 

^=-^°=8^(^+H' ^^=^^'^-^(^+H' ^^=^^^^(^-H' 

These conditions imply that for Rq > 0, both the WEC and SEC are satisfied. DEC implies, on the other hand, from 
(A7) that 

Peff = lEn = 7^ {f - -.Rl] > 0, (36) 



which yields 



3,tl ' 247ri?o 



i?o > and F>^Rl^r< (/^. (37) 




In addition to the energy conditions one can impose the causahty condition (CC) from (A9) 



< 



eff 



< 1, 



(38) 



which is satisfied if F > jRq or r < 



Finally, if one introduces a new parameter (as the equation of state function w) by w = one observes that in 

the range for < r < oo we have 



!<.<-. 



In terms of the physical parameters, if 



then — 1 < a; < 0, and if 



/2g2 



3RI 



< r 



(39) 



(40) 



> r 



(41) 



we have < a; < g. It is clearly seen that the foregoing bounds serve to define possible critical distances where the 
sign of the effective pressure changes sign. This may be interpreted as changing phase for example, from contraction 
to expansion or vice versa in a universe centered by a black hole. We note that scaling the mass and distance by _Ro 
the results will not be affected. For this reason we set Rq = 1. From Eq. (34) we plot in Fig. 2, Cq, (with Rq = 1) 
versus rn and Q. The shaded region for rn < Vc and Cq > 0, which lies below the curve r/, = Vc is the stable region 
outside the black hole. All the rest with Cq < is thcrmodynamically unstable region. Fig. 2 reveals that except for 
a very narrow band of stability islands there is a vast region of instability for Vc at which the effective pressure turns 
sign and continues into opposite pressure, i.e. expansion reverses into contraction or vice versa. 

B. d— dimensions 

1. Thermodynamics 

In higher dimensions one obtains for the entropy and Hawking temperature the following expressions (for n > 2) 



(2ra-l)r: 



„2(n-l) 



2n 



4r(n+i) ' 
4r(n+i) ' u — zn-t 1 



(42) 



Th = 



-1 

„2n-l 



^ ( % (6n - 4n^ + R^rl) rl^^'^Rr' 



2n + l 

8Q~2~ 



47r(2n+l)r2"Jio ^ 

The specific heat also fallows as 



2n- 



- [ ^"-^f"-'^ ) ' +[i?or^-2(2n+l)(n-l)]rr-%" 



d = 2n 



, d = 2n + l 



(43) 
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in which we have used the following abbreviations 

t, = 40" + ,g„ _ + R^^2j rl"-"RS-' (45) 

f, = -4q- - - " V (2„ _ 1) + (_6„ + 4„» + 

= ( '"-"f"-" ) + ^ - 2 (2„ + 1) („ - 1)] -flf^ 

= _8„0-^^ ( "■ + "f ■ - " ) + g2_I) 2 (2,. + 1) (,. - 1)] 

d 

We notice that in odd dimensions from / (i?o) — R§ , i?o can not get negative values for d =odd integer. The detail 
can be seen in Appendix. 



2. The First Law of Thermodynamics 

As it was shown in Ref. Q the first law of thermodynamics in f{R) gravity can be expressed as 

TdS -dE = PdV 

in which E is the Misner-Sharp K§\ energy stored inside the horizon such that 



dE^ — 
2k 



(d-2) (d-3) 



fR + if- RfR) 



Ahdrii, 



(46) 



(47) 



T = ^ is the Hawking temperature, S = ^^^'^ Ir, is the entropy of the black hole P ^ = Tq is the radial pressure 
of matter fields at the horizon and dV — Ahdrh is the change of volume of the black hole at the horizon. In the case 
of constant curvature i.e., R = Rq one gets 



dE= — 
2k 



which implies 



E : 



(d-2)(d-3) d 
rf(d-3) 



4k 



d 

1 

2 



rh 



rh{d+l)Ro d-1 



R§ Ahdrh 



R^Ah- 



(48) 



(49) 



Here we show that the first law of thermodynamics for the metric function (25) is satisfied. Herein P = 



1_ I (d-2)(d-3)Q^ 



and therefore the right hand side reads 



PdV = - — 
47r 



1 /(d-2)(d-3)( 



2 \ 4 



Ahdrh- 



On the other side we have 



TdS-dE = A'^''-^^Rr'drh~^ 
4k rh 2k 



(d-2)(d-3) d 



We combine the latter with (46) and (50) to rewrite the first law as 



2i?n 



1 - 



R§ Ahdrh- 



(50) 



(51) 



4k rh 
or equivalently 



2k 



{d-2){d-S) d 



2Ro 





d 




0-1)1 








47r \ 



(52) 



rh a 



did-2)rt'RS ' 



(d- 2)(i- 3g^ 



2\ 4 



(53) 
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which is the derivative of the metric function at r = . This shows that the first law of thermodynamics by using the 
generahzed form of the entropy for the Misner-Sharp energy is satisfied. To conclude this section of thermodynamics 
we must admit that we don't feel the necessity of addressing the second law. This originates from the fact that we 
are entirely in the static gauge so that the entropy change is assumed trivially satisfied i.e. AS' = 0. 

IV. CONCLUSION 

A relatively simpler class of solutions within / (i?) gravity is the one in which the scalar curvature i? is a constant 
i?o (both i?o > and Rq < 0). We have concentrated on this particular class with the supplementary condition of zero 
energy-momentum trace. The general spherically symmetric spacetime minimally coupled with nonlinear Yang-Mills 
(YM) field is presented in all dimensions {d > 4). The YM field can even be considered in the power-law form in 

d 

which the YM Lagrangian is expressed by L{F) {F°'.F°-)*. Since exact solutions in f {R) gravity with external 
matter sources, are rare, such solutions must be interesting. The equation of state for effective matter is considered 
in the form Pgff = ojp, which is analyzed in Appendix A. The general forms of ui (r) given in (A21) determine u! 
within the ranges of —1 < oj < and < co < ^rj respectively. The fact that a; < — 1 doesn't occur eliminates 
the possibility of ghost matter, leaving us with the YM source and the scalar curvature Rq. In case that the YM 
field vanishes {Q — >■ 0) the only source to remain is the effective cosmological constant Ae/ / = ^'^"2^'"° i which arises 
naturally in / (Rq) gravity. Another interesting result to be drawn from this study is that the effective pressure Pgff 
changes sign before / after a critical distance. Thus, it is not possible to introduce a simple ui =constant, so that the 
pressure preserves its sign in the presence of a physical field (here YM) in the entire spacetime. From cosmological 
considerations the interesting case is when the critical distance lies outside the event horizon. This is depicted in 
the projective plot (Fig. 2) of the heat capacity versus horizon and the charge. Finally it should be added that 
although f{R) = R'^f'^ gravities face viability problems in experimental tests the occurrence of sources may render 
them acceptable in this regard. 
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Appendix A 
Energy conditions 

When a matter field couples to any system, energy conditions must be satisfied for physically acceptable solutions. 



We follow the steps as given in [T^ . 



A. i?o > 



Weak Energy Condition (WEC): 



The WEC states that 



P > 0, 
P + Pi > 0. 



(Al) 



In which p is the energy density and pi are the principal pressure components given by 




I = 




(A2) 



Both conditions are satisfied. So WEC is held. 
Strong Energy Condition (SEC): 
This condition states that 



d-l 



P+EP^ > 0, 



(A3) 



P + Pi > 0. 



The second condition is satisfied but first condition implies that 




(A4) 



or consequently 




(A5) 



By a substitution from (11) for F one finds that for r < rc the condition is satisfied in which 




(A6) 
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Dominant Energy Condition (DEC): 

In accordance with DEC, the effective pressure must not be negative. This amounts to 



1 1 Ro F-^ (rf-2)(rf-l) \ 
= d^lSi ' = (d^2^ 8 ) - ^' 



which for r < it is fulfilled in which 



'(d- 2) (d- 3)( 



{d - 2) (d - 1) Y 4i?j 



Causality Condition (CC): 

In addition to the energy conditions one can impose the causality condition 



This is equivalent to 



which for r <fc ^ satisfied. 

Finally we introduce w = — given by 



(d-2)(d-l) 



which is bounded as 



(d-2) 



d- 1 

It is observed that 



< a; < -rir if r < f 



1 

l<w<0 iffc<r 



(A8) 



0<^ = -^ ] ^<1- (A9) 



^.^-^_(d-2)(d-l)^^ (AlO) 



(All) 



- 1 < w < (A12) 



(A13) 



B. i?o < 

As one may see, presence of Rq in the definition of p and Pi imposes that 7^ 2n + 1 where for n = 2, 3, 4, .... For 

d = 4n we get 

_ -0 _ (F^ ^ 2n-l 

" ^ " 8^n U'^-l^o" ' ' ' 
_ -1 _ |i?o| / 2n- 1 

" " 8^^ + 
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WEC: These expressions reveal that the condition p > and p + > are not satisfied. Similarly the SEC is also 
violated and since the source is exotic we shall not consider it any further here. A case of interest for i?o < is the 
choice = 4n + 2 for n = 1, 2, 3, ... in which 

- "^°-47r(2n+l) 1,^^^ " 2 J ' ^^^^^ 



47r(2n+l) I 2 



WEC: p>0 yields 



2n + l 

r 2 n 

T-TT - - > (A16) 

or 

r < fc (A17) 

where 



2 /n(4n- 1)02 



(A18) 



SEC: The conditions are simply satisfied. 
DEC: This amounts to 



Peff = — ( + !^ + 2nA > 0, (A19) 

"^•^ 4n + 147r(2n + l) I ^ 2 i" ' ^ ^ 



which is also satisfied. 

CC: The causality condition implies 



or equivalently 



which is satisfied for 



where 



0<^= ^ ^ < 1, (A20) 



(4n + 1) 



F 2 _ n 



|i?or+^^<F^ (A21) 



r < rc (A22) 



Here the state function w = becomes 

p 



12 



(4n + l)(^-i 

which is bounded as 



2n + l 

= ^ V. (A24) 



One can show that 



- 1< w < -^—r . (A25) 



-1< w < if < r • ^^^^> 



Figure Captions: 

Fig. 1: The plot of horizon radius rh in 4— dimensions versus mass m for different charges, Q = (Fig. lA) and 
(5 = 1 (Fig. IC). We also plot the heat capacity Cq versus the horizon radius for Q = (Fig. IB) and Q = 1 (Fig. 
ID). Fig. ID, displays in particular the instability caused by the nonzero charge. 

Fig. 2: The 3-dimensional picture of Cq versus rh and Q as projected into the {rh,Q) plane. The shaded region 
with Cq > shows the thermodynamically stable region. From cosmological point of view the region of interest is 
when the critical Tc is outside the event horizon. As shown, below the cruvc = Tc wc obtain stability (dark) regions. 
Above the curve rh = rc, the region is already inside the black hole and no stability is expected. 
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